Solar nanoflares and other smaller energy release 
events as growing drift waves 



J. Vranjes and S. Poedts 
Center for Plasma Astrophysics, and Leuven Mathematical Modeling and 
Computational Science Centre (LMCC), Celestijnenlaan 200B, 3001 Leuven, Belgium. 



Abstract: 

Rapid energy releases (RERs) in the solar corona extend over many orders of mag- 
nitude, the largest (flares) releasing an energy of 10^^ J or more. Other events, with a 
typical energy that is a billion times less, are called nanoflares. A basic difference between 
flares and nanoflares is that flares need a larger magnetic field and thus occur only in ac- 
tive regions, while nanoflares can appear everywhere. The origin of such RERs is usually 
attributed to magnetic reconnection that takes place at altitudes just above the transition 
region. Here we show that nanoflares and smaller similar RERs can be explained within 
the drift wave theory as a natural stage in the kinetic growth of the drift wave. In this 
scenario, a growing mode with a sufficiently large amplitude leads to stochastic heating 
that can provide an energy release of over 10^^ J. 

PACS numbers: 96.60.P-; 96.60.qe; 52.35.Kt 
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I. INTRODUCTION 



The central problems regarding rapid energy releases (RERs) are explaining i) where 
the energy comes from, ii) by what mechanism it is released, and iii) how it is transformed 
into (random and directed) particle energy. A widely (but not universally) used concept 
in the explanation of flares and other RERs, is that of magnetic reconnection which 
implies the presence of electric currents. Yet, the origin, and consequently the source 
of these currents, remains not well understood. Most frequently, shearing photospheric 
motions and/or magnetic flux emergence are used [Ij in the models. The concept has 
been questioned in the past [21 E] because of lack of a measurable magnetic energy and 
configuration change in some events [2 Hj. 

The electric field associated with RER events has been studied extensively in the 
past [SI El El E], with very strong field reported pj of 7 ■ 10^ V/m , that may go [6j up to 
1.3-10^ V/m. As a rule, its presence is associated with the acceleration of lighter particles. 
Apart from the particle acceleration, flares and other smaller energy release events are 
responsible for strong local heating [9] . However, within the current models, it is nontrivial 
to explain how the released energy is dissipated so fast in the highly conductive corona, 
characterized by a Lundquist number around 10^^. The existing continuum approximation 
(fluid) or magnetohydrodynamics (MHD) models, are not very promising because it is 
clear that the actual heating takes place at length scales much smaller than those on 
which the (macroscopic) MHD model is justified. We stress also that the possibility can 
not be excluded that, at least in some cases, the reconnection can appear as a consequence, 
rather than a cause of a RER. Another problematic aspect regarding the role of RERs 
in the heating, is the distribution of magnetic null points. According to Ref. [TO] only 
2% of them are located in the corona and 54% in the photosphere. This is opposite to 
what would be required for the mechanism that is presently believed to heat the corona. 
The heating by waves rather than by reconnection is also supported by the diagnostic of 
active regions presented in Ref. 

Here we present an alternative model of RERs and the consequent plasma heating, 
based on the kinetic theory of drift waves, driven by density gradients that are omnipresent 
in the solar corona. The model is able to provide answers to all the fundamental ques- 
tions raised above. The required density gradients are visible in the magnetic loops 
spread throughout the solar atmosphere. They are closely connected to the continuous 
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restructuring of the magnetic field (due to the frozen- in conditions), and also to the ob- 
served inflow of the plasma along the magnetic loops [12j. Measurements by Voyager 1 
and Voyager 2 show [T3l [T4] that the transverse size of some of these highly elongated 
structures at the Sun can be below 1 km. Hinode observations |15] only confirm that 
the solar atmosphere is a highly structured and inhomogeneous system. Moreover, a 
three-dimensional analysis \16\ of the coronal loops reveals short-scale density irregular- 
ities within each loop separately. Presently, the observable characteristic dimensions of 
density irregularities are limited by the available resolution of the instruments (at best 
about 0.5 arcsecond). However, even extremely short, meter-size scales can not be ex- 
cluded. This can be seen by calculating the perpendicular (to the magnetic field) diffusion 
coefficient for a corona environment Dj_ j ^ nTjVj / (nijVt^) oc m^j'^, j = i, e. The diffusion 
velocity in the direction of the given density gradient is [T7] D^ jVno/nQ. Taking the 
inhomogeneity scale-length L„ = [{dno/dx)/nQ]^^ = 10, 10^, 10^ m, where x-denotes the 
direction perpendicular to the magnetic field vector, we obtain the ion diffusion veloci- 
ties, respectively, 10^'^, 10^^, 10^^ m/s only. Therefore, even very short density inhomo- 
geneities can last long enough to support relatively high frequency drift instabilities. In 
dealing with the drift wave, we may thus operate with the density inhomogeneity scale 
lengths that can have any value from meter size up to thousands of kilometers in the case 
of coronal plumes. The role of the drift wave in RERs has been overlooked so far in the 
literature, probably due to the fact that it simply does not exist in the widely used MHD 
model. 

The purpose of this work is to show that the free energy stored in these density 
gradients may drive drift waves, and these may further release energy on a massive scale. 
The dissipation of the drift waves is easy to explain in our self-consistent kinetic model 
that works on the (very small) length scales at which the actual dissipation takes place. 
Actually, two mechanisms of energy exchange and heating will be shown to take place 
simultaneously, one due to the Landau effect and another one, stochastic heating which 
acts exclusively in the perpendicular direction. To prove this we use only established basic 
theory, verified experimentally in laboratory plasmas. 

II. MODEL AND RESULTS 

The drift wave properties within well known limits are described |T8] by the frequency 
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and the growth rate, respectively 
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Here, Ao(&i) = /o(&i) exp(-6i), 6^ = fc^pf, = t-Ti/f^i, A^, = Vn/ujpi, uj*i = -u^eTi/Te, 
= kyV^ei '^*e = —{v'^e/ ^e)(^z X V^Tio/no, Jq is the modified Bessel function of the first 



In Eq. (j2]) the positive second term in the square bracket describes the damping on 
ions. The necessary condition for the instabihty Ur < oj^^e-, from the first term is, as a 
rule, easily satisfied. Note that Eq. ([1]) reveals the presence of the energy source already 
in the real part of the frequency Ur oc V^rio, while details of its growth due to the same 
source are described by Eq. ([2]). 

In Fig. 1 we give the wave growth-rate as a function of two parameters, the perpen- 
dicular wavelength and the density scale-length L„. The chosen parameters correspond 
to the inner corona, i.e., Bq = 5 ■ 10"^ T, no = 10^^ m-^ T, = Te = To = 7 ■ 10^ K. The 
graph does not change drastically by varying these parameters for ±1 order of magni- 
tude. The wavelength parallel to the magnetic field vector is taken a.s Xz = s ■ 10^ m. We 
have introduced a parameter s that we can vary within the range 10^^ — 10^ in order to 
demonstrate that the mode behavior remains unchanged at various density scale lengths 
Ln, provided that we keep the ratio Xz/Ln constant. It can easily be shown that in this 
case the ratio 'j/oOr also remains constant [19l|20], while both quantities are shifted either 
to lower or higher values. As example, the mode frequency for s = 1, L„ = 100 m, and 
Xy = 1 m is uJr = 231 Hz, and it changes as 231/s when the other parameters are fixed. 
The growth rate may easily become larger than Ur for short Xy,Ln- 

The growth of the wave results in a stochastic heating mechanism that implies single 
particle interaction with the wave. The process has been experimentally verified [211 122] . 
In the drift wave, the ions move in the perpendicular direction to large distances and they 
feel the time- varying field of the wave due to the polarization drift Vp = {dE / dt) / (QiBo) , 
and as a result their motion becomes stochastic. This stochastic heating is anisotropic, 
acting mainly in the direction normal to the magnetic field Bq. The perpendicular heating 



kind and of the order 0, Bq = BoCz, and V_l = Cxd/dx. 
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Figure 1: The growth rate ([2]) normahzed to the wave frequency Ur in terms of the 
perpendicular wavelength and the density scale-length L„. 

in the experiment was larger by about a factor 3 compared to the parallel one [211 122] , 
and it is exceptionally fast. It the same time, in view of the mass difference and the 
given physical picture, this heating scenario predominantly acts on the ions, with heav- 
ier ions more efficiently heated than lighter ones. For drift wave electrostatic potential 
perturbation of the form (l){x)sm{kyy + kzZ — ut), \ky\ ^ \kz\, one finds the ion po- 
larization drift velocity following the procedure from Ref. [23]. For large enough par- 
ticle displacements (large wave amplitudes) Vp = ey[ky(f)/{QiBo)]{d/dt) sm{kyy — ut) = 
ey[ky(f) / {QiBo)]{kydy / dt — uj) cos{kyy — ut), where dy/dt = Vp. This yields 

ujky(j) cos{kyy - cut) 
^ ^ ^liBo 1 — a cos{kyy — cut) ' 

a = kyp'^e(f)/{K,Ti). 

It has been shown [22] that the stochastic heating takes place for a large enough wave 
amplitude, more precisely for a > 1. This condition implies that the ion displacement due 
to the polarization drift is comparable to the perpendicular wavelength. This is because 
Vp = X dv±/dt)/Qi, and v± is the leading order E x 5-drift, so that Vp ~ au/ky and 
the perpendicular displacement due to the polarization drift is 5 = Vp/u = a/ky. Another 
important feature is that Vp ~ ky, hence the stochastic heating is due to the electrostatic 
property of the wave. 
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Figure 2: Volumetric energy released during stochastic heating by the drift wave in terms 
of the perpendicular wavelength. 

According to [2n [22] . the maximum achieved bulk ion velocity, proportional to the 
wave amplitude, is given by 

~ [klple<P/{KT,) + 1.9]a/fc,. (4) 

As explained in Ref. [22], the factor 1.9 appears after making Pincare plots of particle 
trajectories for different values of a. Slightly different values of this factor, i.e., 1.5 and 
2.3, are reported in Refs. [2HI25], respectively. 

Eq. (jl]) implies an effective increase in the temperature Te/ = rriiv'^/ (?>k), while the 
volumetric energy increase of the stochastically heated particles is = riQUiiv'^ [J/m^], 
and the energy release rate is Vm = ^mj^g [J/ (m^s)], where Tg = I/7. 

Equation (jl]) reveals a minimum of Vm at Xy = 27rpj[e0/(fi;Tj)]^/^. This is seen also in 
Fig. 2. Here, we fix = 10^ V and use the same parameters as in Fig. 1. 

Since the results are valid for any s, the growth time oc s may have rather different 
values. Assuming a small accidental initial perturbation e(f)o/{i^Ti) = 0.01 the growth 
time till it gets the value used above is Tg = log(0/0o)/7- Taking s = 1 and s = 10^, at 
the minimum from Fig. 2 we obtain Tg ~ 0.2 s, and Tg ~ 3.3 min, respectively. The same 
dependence on s holds for the total plasma volume Vp = LxLyLz involved in RERs, and 
consequently for the total amount of the released energy. We may take = Ln, and for 
Lz take one wavelength Xz only while for Ly we may take a layer of around 10 km. At 
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the minimum A^, from Fig. 2 and for s = 10^, this yields I^mVp = 1-5 ■ 10^^ J. Clearly, 
this value can easily become considerably larger (e.g., for a larger volume, then for A^^ 
not taken at minimum, and also increasing the density for one order of magnitude will 
increase the energy for one order too). 

Calculating dTf,f{mi)/dmi it can be shown that the heating increases with ion mass if 
kypf[e(j)/ (uTi)] > 1.9. For the parameters of interest here, it turns out that this condition is 
always satisfied. We stress that numerous indications and observations in the solar corona 
[26t [271 [28] confirm these peculiar features of a stronger heating of heavier particles, i.e., 
Tg < Th < Thc, and a stronger heating in the perpendicular direction as compared to the 
heating along the magnetic field vector. 

The drift mode presented above implies a time-varying electric field. Its parallel 
component |V^0| for A^ = s ■ 10 km (in units of the Dreicer runaway electric field) is 
4 and 0.004 for s = 1 and 10^, respectively. In the first case, the bulk plasma species 
(primarily electrons) can be accelerated by the wave in the ^-direction. In the latter, this 
holds only for electrons from the tail in the distribution function. Hence, in any case there 
will be an increase in the directed electron energy. The escaping electrons imply that a 
smaller amount of them are available to shield ion perturbations and the mode should be 
even more growing. The electron parallel velocity in a time-varying parallel wave-electric 
field Eq cos{kzZ — Urt) is 

eEz 

v{t) =Vo {sm[k^Zo + {k^vo - uJr)t] 

-sm{k^Zo)}. (5) 

Here Vq and zq are the starting electron velocity and position in the parallel direction, 
respectively. A strong acceleration will take place, in particular for particles close to 
resonance vq = oJr/kz = O.lvre- The non-relativistic energy radiated by an electron decel- 
erated in the time- varying wave electric field is given by AS/At = {Av / At)"^ / {GTceoc^) ■ 
Here, At ~ Tosc/4, where Tqsc = '^'^/{kzVo — oJr)- It is seen that for particles close to reso- 
nance, the oscillation period Tq^c becomes very large and the same holds for the particle 
velocity. So the previously described process of stochastic heating will be accompanied 
by a large directed acceleration and by radiation as well. 

For plasma (3 above the electron/ion mass ratio the drift wave couples with the Alfven 
wave. This couphng is given by [29l [30] a;^ — uj'^{uj^e + + uj[u)^e^^i — k^cl/ (1 + kypj) — 
k^k'^cKpf + pD] + klcluj^e/i'^ + k^pf) = 0, = Bl / {poUomi) . It describes the drift 
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wave and two Alfven waves. The electromagnetic effects in RERs are frequently (but not 
always) observed. However, even in the presence of coupling, the previous analysis related 
to the energy release will not change considerably. This can be checked by solving the 
above dispersion equation for the coupled modes and comparing to previous results. As 
example, for Xy = 3 m, s = 10^ and other parameters unchanged, the drift wave frequency 
from Eq. ([T]) is 0.21 Hz, while from the coupled mode it is 0.25 Hz. 

III. SUMMARY 

The model and results presented here, based on the drift wave theory, yield an alter- 
native description of some rapid energy releases in the solar corona. Because of a high 
temperature and a low number density, collisions are rare in coronal magnetic structures, 
hence the kinetic description of the drift wave is the most appropriate. In such a descrip- 
tion the drift wave is strongly growing and its growth-rate appears as a purely kinetic 
effect. 

Using parameters typical for inner corona we have shown that such a growing drift 
wave can lead to stochastic heating of ions,^'' with the heating rate so large that can be 
used for the description of nano-fiares and smaller energy releases. We have also shown 
that there should be an acceleration of plasma particles, primarily electrons, by growing 
drift waves in solar corona. Such an acceleration and heating of coronal plasma has 
been discussed in many studies in the literature in the past,^^~^^ yet not within the drift 
wave theory. Our analysis comprises the electrostatic limit, which may be appropriate at 
least in some cases because observations in the past have shown that not all rapid energy 
releases are associated with a measurable change in the magnetic field topology. However, 
using the well known theory the magnetic effects can be included and they would lead to 
the coupling between the drift and Alfven waves, and this limit is also briefly discussed. 

The model implies the density gradients in the background plasma, that are expected 
within coronal magnetic structures. Since the growth of the mode is on the account of 
the energy stored in this gradient, it is expected that the background is simultaneously 
changed in the process of the instability. Simulation in the past have shown this.^^ As 
a result the energy release should be below the value obtained in the present study. A 
detailed description of this feedback effect can be performed only numerically, yet this is 
beyond the scope of the present work. 
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